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Navier-Stokes turbulence subject to solid-body rotation is studied by high-resolution 
direct numerical simulations (DNS) of freely decaying and stationary flows. Setups char- 
acterized by different Rossby numbers are considered. In agreement with experimental 
results strong rotation is found to lead to anisotropy of the direct nonlinear energy flux, 
which is attenuated primarily in the direction of the rotation axis. In decaying turbu- 
lence the evolution of kinetic energy follows an approximate power law with a distinct 
dependence of the decay exponent on the rotation frequency. A simple phenomenologi- 
cal relation between exponent and rotation rate reproduces this observation. Stationary 
turbulence driven by large-scale forcing exhibits fcj^-scaling in the rotation-dominated in- 
ertial range of the one-dimensional energy spectrum taken perpendicularly to the rotation 
axis. The self-similar scaling is shown to be the cumulative result of individual spectral 
contributions which, for low rotation rate, display fcj^-scaling near the fcy — plane. A 
phenomenology which incorporates the modification of the energy cascade by rotation is 
proposed. In the observed regime the nonlinear turbulent interactions are strongly influ- 
enced by rotation but not suppressed. Longitudinal two-point velocity structure functions 
taken perpendicularly to the axis of rotation indicate weak intermittency of the /cy = 
(2D) component of the flow while the intermittent scaling of fcy 7^ (3D) fluctuations is 
well captured by a modified She-Leveque intermittency model which yields the expression 
Cp — P/6 + 2(1 — (2/3)^/^) for the structure function scaling exponents. 



1. Introduction 



Hydrodynamic turbulence subject to rotation is an ubiquitous problem in fiuid me- 
chanics. The understanding of its detailed properties is crucial for engineering problems 
such as the design of turbomachinery and for the understanding of atmospheric and 
oceanic flows influenced by the earth's rotation. These applications have motivated an 
extensive research on the macroscopic and spectral properties of rotating turbulence 
though a comprehensive and fully consistent physical picture is still missing. 

Early experiments show the tendency of rotating turbulence to asymptotically two- 
dimensionalize in planes perpendicular to the rotation axis il, where f2 — \^l\ is the 
rotation frequency. Signatures of this behaviour are, e.g., an increased ratio of velocity 
correlation lengths along and perpendicular to the axis of rotation, (see, for example, Ib- 



betson & Tritton 1975 ) , visible elongation of spatial structures along n, (e.g. in Hopfinger 
et aL||1982 ), or anisotropy of characteristic length- and time-scales (e.g. Wigeland||1978 



Jacquin et al. 19901. In addition, slower decay of kinetic energy is found in the works 



of Wigeland ( 1978 1 and Jacquin et al. ( 1990 1 as compared to non-rotating turbulence. 
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More recently, the experiments of Morize & Moisy ( 2006 ) show the significant influence 



of confinement on the decay of rotating turbulence. Baroud et al.\ (2002, 2003) consider 
two-point increment statistics up to order p = 10 yielding inertial-range scaling exponents 
Cp = The work of Morize et al. (2005') focuses on the energy spectrum E{k) which 
exhibits a Rossby number dependent scaling exponent continuously ranging between -5/3 
(slow rotation) and about -2.3 (fast rotation). Both experiments feature turbulent flows 
with energy injection at intermediate scales and thus allow an inverse cascade of kinetic 
energy to develop in the quasi-two-dimensional flow expected for strong rotation. 

The first direct numerical simulations of rotating turbulence conducted by |Bardina| 
et al. (19851, Mansour et al. (19921 and Hossain (19941 also suggest reduction of the en- 



ergy fiux and two-dimensionalization. They suffer, however, from small spatial resolution 
and correspondingly low Reynolds numbers and do not yield conclusive results on the 
spectral properties of the fiow. Simulations at higher Reynolds number by |Godeferd k,\ 



LoUini (1999| and Morinishi et al. (j2001j provide more evidence for the experimentally 



observed behaviour of the correlation lengths and the decay properties of the kinetic tur- 
bulent energy, respectively. For the case of driven turbulence Yeung & Zhou ( 1998 1 find 
fc~^-scaling of the energy spectrum in a setup with large scale forcing. On the contrary, 
the simulations of [Smith fc Waleffe| ( [T999| , |Chen et a?.| ( |2005| ) and |Mininni et a?.| ( [2009 l 
explore systems with forcing at intermediate scales. All three works identify a quasi- 
two-dimensional state in the case of strong rotation accompanied by an inverse energy 
cascade for wavenumbers k < fcforcing- In this wavenumber range Smith & Waleffe ( 1999 1 
report to find fc~^-scaling of the energy spectrum while a A:~^-behaviour for k > fcforcing 
is observed by [Smith fc Waleffe, (1999) and Mininni et a?.[(|2009[) 



Large-Eddy simulations like those of Squires et al. (1993), Bartello et al. (19941 and 



Yang & Domaradzki ( 2004 1 are a useful tool to study the effects of rotation on char- 



acteristic properties of rotating turbulence such as the slowing down of energy decay, 
the large-scale structure of the flow and its tendency to become two-dimensional. The 
approach allows to attain Reynolds numbers beyond the reach of DNS. However, the 
results of LES crucially depend on the applied subgrid model which requires careful 
adjustments and gauging by comparison with high-resolution DNS or experimental mea- 
surements. For the detailed numerical investigation of small-scale properties of turbulence 
like anisotropic inertial-range scaling of the energy spectrum and higher-order structure 
functions, therefore, the DNS approach is chosen in this work. 

Different theoretical approaches exist with regard to rotating turbulence. Shell- models, 
which are strongly simplified representations of nonlinear turbulent dynamics, allow for 
the straightforward inclusion of rotation effects, e.g. via modes governed by a time- 
correlated stochastic process as in jHattori et al. (2004). The simplicity of those models. 



however, restricts their predictive capabilities to a continuous change of spectral scaling 
between Kolmogorov fc~^/'^-scaling (slow rotation) and fc~^-behaviour for fast rotation. 
Zhou (1995) and Mahalov & Zhou (1996) recover fc~^-scaling of the energy spectrum in 



the context of quasi-normal closure theory using dimensional arguments first employed 



by Kraichnan (1965). Canuto fc Dubovikov (1997) arrive at similar results via a formal 



treatment of the spectral energy flux using helical-mode decomposition. The former work 
is based on the assumption that tq = is the dominant time scale of nonlinear relax- 
ation while the latter study is analytically showing the same covering cascade dynamics 
for rotation rates up to tq ~ tnl i/vi. Here, tnl is the nonlinear eddy-turnover or 
relaxation time defined with the RMS velocity ve at scale i (see below). On the con- 
trary, weak turbulence theory (see, e.g., Cambon et al. 2004 Galtier 2003 [Bellet et al. 



2006 1 yields different anisotropic scaling relations but is only strictly applicable in the 



asymptotic limit tq <C tnl- 
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This article is motivated by the lack of a simple dynamical picture of spectral en- 
ergy transfer which could also help to improve the understanding of turbulent energy 
decay. The proposed phenomenology is backed-up by high-resolution direct numerical 
simulations of incompressible rotating homogeneous turbulence in free decay and of ro- 
tating turbulence subject to large-scale driving. The kinetic energy is found to exhibit 
approximate power-law decay with an exponent that decreases with increasing rotation 
frequency. This behaviour is in agreement with the observed attenuation of nonlinear 
spectral transfer under the influence of rotation and is captured by a simple model 
based on the cascade phenomenology. In the case of forced turbulence, analysis of one- 
dimensional spectral data indicates anisotropic energy flux mainly perpendicular to il. 
Inertial-range scaling of the perpendicular energy spectrum ~ fcj^ is observed and re- 
produced by the dynamical cascade phenomenology. While the parallel spectrum does 
not show clear scaling, the perpendicular spectra at fixed A:|| display fc~'^-scaling in the 
vicinity of = for low rotation rate. For growing il, energy accumulates at fc|| = in 
agreement with the observed quasi-two-dimensionalization of the system. The intermit- 
tency measured via velocity structure functions perpendicular to ft is very weak for the 
2D component of the flow (/cy =0). In contrast, the intermittent scaling signature of the 
3D component (fc|| ^ 0) is in agreement with a modified She-Leveque model and exhibits 
a weak trend towards Gaussianity for high rotation rate. The remainder of this paper is 
organized in the following way: In section 2 we introduce the equations of rotating fiuid 
flow and present the numerical methods that have been used. Sections 3 and 4 comprise 
the results on spectral and macroscopic properties of the flow and a summary is given in 
section 5. 



2. Basic equations and numerical methods 

Incompressible fluid flow subject to solid-body rotation is usually described in a frame 
of reference rotating about a fixed axis ft = fle^ with frequency CI using the dimensionless 
Navier-Stokes equations including the Coriolis force. 



dtv + V -Vv ^ -~VP + fiAv + 2flv X e^, 
V • t; 0, 



(2.1) 
(2.2) 



where the centrifugal force has been incorporated in the generalized pressure P, or the 
numerically more favourable vorticity formulation 



dfUj = Vx {vxuj + 2flvx Bz) + ^Aw, 
\7 ■v = 



(2.3) 
(2.4) 



with the vorticity u) — V x v (see, e.g., Greenspan 19681. The dimensionless kinematic 



viscosity fi and rotation rate are both assumed to be constant. 



Equations (2.3 1 and (2.4) are integrated numerically on a cubic box extending 27r in 



each dimension with triply periodic boundary conditions. The integration is performed 
in a pseudospectral Fourier representation by applying an explicit leapfrog scheme. The 
diffusive term is incorporated using an integrating factor technique (see, e.g., Meneguzzi 
fc: Pouquet|1989 1. The time step which is required for numerical stability by the Courant- 
Friedrichs-Lewy criterion ensures the temporal resolution of all inertial-wave oscillations 
present in the system. The aliasing error introduced by the pseudospectral approach is 
reduced by spherical mode truncation. The quasi-stationary state of the driven systems 
is sustained by a forcing which freezes all modes with k ^ kf = 2. The applied resolution 
is 512'^ collocation points for forced and 256"^ for decaying turbulence. 
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simulation grid forcing fj, tsim Ro Re 

I 512^ yes 5 4 x 10"* 45 5.3 x 10"^ 4000 

II 512^ yes 50 4 x 10~* 30 1.3 x 10"^ 2300 

III 256^ no 0-5 1 x 10"^ 36 oo - 2 x 10~^ 300-1100 

Table 1. Parameters for the different simulation runs: rotation rate Q,, viscosity /x, total 
duration tsim, Rossby number Ro, and Reynolds number Re 



The isotropic initial state of the forced simulations is a smooth vorticity field with 
random phases and an energy spectrum E{k) ~ exp (— fc^/fco) , fco = 4. This configuration 
is integrated forward in time without forcing until the maximum of enstrophy is reached. 
After this period of time, which corresponds to one large-eddy-turnover time (LET) and 
to three units of time in the notation of this paper, the energy spectrum is fully developed 
in its spectral extent and the forcing is switched on. The spectral distribution of energy 
evolves self-consistently towards a state with a low-order power-law at largest scales and a 
Kolmogorov-type inertial range. As soon as total energy E = Jy dVv'^/2 and dissipation 
e = jjL jy dVuP' are statistically stationary with ~ 1 (more precisely: 0.87 for = 5 
and 1.3 for = 50) and e ~ 0.05, is set to a finite value. For decaying turbulence 
simulations a similar random initial field is generated and left to evolve freely for one 
LET before rotation is switched on. After the onset of rotation a statistically stationary 
state of turbulence is maintained in driven simulations for 10 or 15 large-eddy turnover 
times depending on the rotation rate. Decaying runs extend over approximately 12 LETs. 
The dimensionless kinematic viscosity is constant in all simulations. 

The characteristic length Lq and velocity Vq , necessary for the calculation of the macro- 
scopic Rossby number, Ro = Vb/(2riLo), and Reynolds number. Re = LoVb//i, can only 
be determined a posteriori in homogeneous turbulence. Therefore, both quantities are es- 
timated in the dimensionless system using E, e, and as Lq ~ E/inef/'^ and Vb-i?^ 
Hence Rossby and Reynolds numbers given in this paper are defined as Ro = y/e/JmE) 
and Re = E^ / {^e) / fi, respectively. Specific parameters of the performed simulations 
are summarized in table [T] 



3. Structure of the flow 

An immediate impression of the structure of rotating turbulence is provided by visual- 
izations of the fiow field given in figure [T| The system shows a distinct structuring along 
the rotation axis which is especially pronounced for simulation II and indicates a strong 
velocity correlation along this direction. The observed trend to a quasi-two-dimensional 
state with faster rotation is in agreement with previous experimental findings (see above). 



It is the consequence of a nonlinear analogue of the Taylor-Proudman theorem (see Ma- 
halov fc Zhou|p96l |Chen et a?.|2005[ ). For more mathematically oriented works on, e.g., 



the dynamics of two- (fcy = 0) and three-dimensional (fcy ^ 0) flow components as well 
as the issues of their non-vanishing mutual interaction in the limit oo see |Babin| 



et at. (20001 and references therein. 



The clearest signature of the turbulent energy cascade is found in energy spectra 
and energy flux functions. Due to the above-mentioned anisotropy of the turbulence, 
one dimensional spectral quantities are considered in directions parallel and perpendic- 
ular to the rotation axis. The axis-parallel nonlinear energy flux over wavenumber fcy 

is given by T^^^ — /g*^" d/cy J dfc' J dfc_L {ioj* ■ (k x [v x u:]^ + c.c.^ /k^, with [•] denoting 
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Figure 1. Visualization of the magnitude of the velocity field \v{r,t)\ for simulations I (left) 
and II (right) in the statistically stationary regime of the flow. Brighter colours correspond to 
higher values of |?;| and darker colours to lower values. The rotation axis points in the vertical 
direction. The pictures show the flow fleld of the whole 512'^ simulation domain. 



Fourier transformation, k± _L fJ, fc|| || ft, and k' _L fe_L,fe||. The axis-perpendicular flux 

function T^^ is defined analogously using the integration J^^ dfc^ / dk' J dfey with the 
fc_L-component running along an arbitrary fixed direction in the axis-perpendicular plane. 
The fluxes are normalized with the total energy dissipation rate e — ^ J dVcu^. 

The parallel and perpendicular energy fluxes for = 0,5, 50 are negative at all 
wavenumbers indicating direct energy transfer towards small scales ( jMiiller fc Thiele 
20071. The two-dimensionalization visible in figure [l] is not in contradiction with the 
observed direct fluxes of energy which are a mere consequence of the large scale driving 
of the flow. The fluxes weaken with increasing fi, a trend that is much stronger in Tk^^ 
than in T^^ . The increasing anisotropy of the energy cascade with growing f2 reflects 
the nonlinear nature of the dynamical process that is at least partially responsible for 
the observed two-dimensionalization of the flow. The classical picture of energy transfer 
through eddy breakup relates a stronger eddy scrambling perpendicular to the rotation 
axis than in the parallel direction along which coherent structures remain comparatively 
intact. 

The perpendicular energy spectra E{k±) = J d/cy / dk' |ufc|^ /2 shown as solid lines in 
figure |2] display inertial range scaling for 4 < < 20 (I) and 5 < fc_L < 16 (II). The 
shortening of the scaling range with increasing rotation rate is due to an extension of 
the dissipation region towards smaller k (see below). In simulation II a bump persists 
around fc/, where the forcing region descends into the freely evolving range of scales. It is 
a result of the simple forcing scheme and is tolerable at the chosen numerical resolution 
as it is not significantly perturbing the flow beyond k ~ 5. 

In spite of the shortened inertial range in simulation II both cases clearly exhi bit the 
scaling E(k±) ~ fc 7^. These results a gree with the direct numerical simulations of |Yeung" 
fc Z hou|(|1998[) and Smith & Waleffe|(|1999|)(for k > kf) as well as shell-model calculations 



(see Hattori et aL\\2004: \. Recent experimental results on decaying turbulence subject to 
rotation (see Morize et aL||200^ l yield an energy scaling exponent w —2.5 for the micro 
Rossby number of Roi^ = {lj'^)'-^'^ /{2fl) ~ 0.08 (corresponding to our case II) and an 
exponent of « —1.7 for Ro^ ~ 0.7 (corresponding to I). The experimental configuration 
with an initial excitation of the flow followed by decay under rotation is however not 



6 



M. Thiele, W.-C. Miiller 




Figure 2. Perpendicular two-dimensional spectra _B(fc||,fcj_) for n = 5 (left) and f2 = 50 (right) 
at fixed fc|| and compensated with k\. The spectrum for A;|| = (dotted) differs markedly from 
those for fcy =1 — 6 (dashed lines from top down) and higher. Three-dimensional perpendicular 
spectra are obtained by fcy -integration from E{k\\,k±_) (solid). The dash-dotted line indicates 
fej^-scaling. 



directly comparable to the simulations described here, which are characterized by a 
continuous large-scale forcing. We note in passing that these Rossby numbers place the 
present simulations in the "intermediate" Rossby number regime of |Bourouiba fc Bartello| 
([2007V In the experiment of [Baroud et al.\ ( 2002| the same scaling ~ A: ^ is observed 
although here forcing at the small scales produces an inverse energy cascade. 

Several theoretical models arrive at the scaling shown in Fi g. [2] In the context of 
quasi-normal closure theories Zhou ( 1995 1 and Mahalov & Zhou ( 1996 ) have conducted a 
dimensional analysis of the energy flux terms, assuming r* ~ To ~TF^ for the relaxation 
timescale of nonlinear interactions r*, whereas Canuto fc Dubovikov| ( 1997 1 have analyzed 
the energy flux formulated in helical mode decomposition, both arriving at the same 
scaling result. On the contrary, weak turbulence theory yields E{k) ^ k^^ (see Cambon 
et aL||2004[ [Bellet et a?.|2006[ ) in the asympto tic quasi-normal Markovian approximation 



valid for /cy ^ and E{k) ~ k~^ (see 'Galtier 2003). The former result is regarded as the 
consequence of a strongly anisotropic spectrsil energy distribution with respect to Q, while 
the latter explicitly assumes spectral isotropy in combination with a phenomenological 
nonlinear transfer time based on the wave-kinetic equation. Generally, however, wave 
turbulence theory predicts anisotropic scalings in the context of rotating turbulence. 

It is noteworthy that the axis-parallel spectra shown in flgure [3] do not show clear 
scaling. This indicates that the perpendicular domain of wavenumbers provides the dom- 
inant contribution to the angle- averaged scaling E{k) ~ fc~^ found in some of the above- 
mentioned numerical and experimental studies. 

The contributions to the perpendicular spectrum for flxed fcy, shown in flgure [2] exhibit 
for = 5 fcj^-scaling in the vicinity of /cy = 0. Note that E{kj_) ^ J dk\\E{k\\, kj_) 7^ 
-E(fc|| — const, k±). The scaling range, however, shrinks from the low k side with growing 
fcj|. For = 50 the spectral contributions for flxed /cy show no clear scaling. Evidently, 
the observed fc7^ power-law is due to spectral contributions from all wavenumbers fcii. 



This observation is in contradiction to weak turbulence simulations reported in Bellet 



et al. 



k 



(2006) which give kj^ -scaling in directions perpendicular to Q, and suggest that 
spectra are generated by spherical averaging. Furthermore and especially for f2 = 50, 



the energy is mainly concentrated around the plane fcii = as expected from previous 



Rotating homogeneous turbulence 



7 




M 10 - 



10 



100 



Figure 3. One-dimensional energy spectra parallel to the rotation axis compensated with fcy 

for n = 5 (sohd), and H = 50 (dotted). 



studies (see, e.g., Cambon fc Jacquin||1989 Cambon et aZ.|[l997 ) and in agreement with 
the observed nonlinear two-dimensionalization. 

The observed axis-perpendicular spectral scaling can be described by a simple phe- 
nomenology of rotating turbulence proposed in Miiller & Thiele (20071 which tries to 
capture the effect of inertial oscillations on convective fluid motion. The phenomenol- 
ogy applies in the inertial-range under the additional constraint fcn > fc > min(fcd, fcj )j 
i.e. when nonlinear energy transfer is dominated by Coriolis force effects. This spec- 
tral interval is bounded by the rotation wavenumber fcji — (f^'^/e)^^^ whose inverse, also 
called Rossby d eformation radius, is analogous to the geophysical Ozmidov scale (see e.g. 

(£/m3-)1/4^ 



Ozmidov 1992), the classical Kolmogorov dissipation wavenumbe r kd 
its rotation-modified counterpart fcj — (e/fi)^/^/^ (see, e.g 



also 



fc Duboviko^[T997 l. 



Zeman 



1994 



Canuto 



For turbulent fluctuations perpendicular to ft on spatial scale ^ kj^ the influ- 
ence of rotation on the inertial-range energy transfer in planes perpendicular to n is 
phenomenologically described similarly to the Iroshnikov-Kraichnan picture of magne- 
tohydrodynamic turbulence (llrosh nikov|| 1964 Kraichnaii||1965 1 . Taking into account the 
random displacements that a fluid particle experiences during a nonlinear interaction 
of turbulent eddies under the influence of inertial oscillations leads to a prolonged en- 
ergy transfer time rtr ~ t^l /'^f2 "^ith tnl^ ~ ^/v^ as compared to the non-rotating 
case. This yields by standard dimensional arguments the non-intermittent scaling re- 
lation ^ (O e)^/^^^/^ corresponding to the observed scaling of the energy spectrum 
~ fcj^ (see also |Zhou||l995| |Canuto fc Dubovikw||l997| . 

The numerical values of the characteristic wave numbers are for simulation (1) : fco — 50, 
fcj ~ 250, kd ^ 167, and for simulation (11): fcn ~ 1581, fcj ~ 79, kd ~ 167. It should be 
noted that, according to their definition, for increasing fi, fcn grows while fcj decreases. 
This is the reason for the smaller inertial range in system 1 and could also explain why 
in earlier works with lower resolution the scaling exponent of the energy spectrum for 
high rotation rates has been difficult to pin down. 

The two-dimensionalization evident in figure [T] due to the large-scale driving, does not 
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lead to an inverse energy cascade but still leaves some spectral traces that can be un- 
derstood by the following ordering considerations: it is proposed that for fixed the are 
three different regimes of rotating turbulence: a) turbulence with negligible rotation ef- 
fects at fc 3> A:si, b) turbulence strongly dominated by the Coriolis force at ^ fcj^ which 
is asymptotically described by weak-turbulence theory with the associated purely non- 
linear two-dimensionalization via anisotropic nonlinear energy transfer, and c) rotating 
turbulence in the transient regime at fc < fcfj observed in the present simulations where 
quasi-two-dimensionalization occurs due to the "nonlinear Taylor-Proudman theorem" 
(see above). This effect entails an observable energetic separation of 2D (fc|| = 0) and 3D 
^ 0) components of the velocity field which exhibit an energetically weaker coupling 
at high rotation rate (see, e.g., Bourouiba fc Bartello^2007, and references therein). It can, 
however, be shown rigorously (see e.g. |Babin et aZ.||2000 l that the dynamical coupling of 
2D and 3D fiuctuations remains finite for all values of f2. In regime a) turbulence is based 
on nonlinear eddy interaction, regime b) represents turbulence governed by (weak) wave 
interaction, while nonlinear dynamics underlying regime c) correspond to wave-modified 
eddy interaction. In regime c) the 2D fluctuations should exhibit dynamics different from 
the 3D fluctuations as indicated by the nonlinear Taylor-Proudman theorem. A possible 
scenario that would be in agreement with the observed scaling laws consists of a direct 
cascade of 3D energy as described by the proposed phenomenology in combination with 
an inverse cascade of 2D energy (towards wavenumbers smaller than the forcing range) or 
a direct enstrophy cascade (towards wavenumbers larger than the forcing range) . In this 
picture the inertial-range scaling of 2D energy could be interpreted as the signature 
of a direct 2D enstrophy cascade. 

Regime c), in fact, describes a state of continuous transition from Ro — > oo (negligible 
rotation) to Ro — > (weak inertial- wave-turbulence). This transitional state is applying 
to most realistic rotating systems. The dynamical phenomenology by [Miiller k, Thiele 



(20071 could be a reasonable starting point for extending the theoretical description 
of rotating turbulence to include the growing influence of inertial waves on turbulence 
dynamics. 

To obtain the scaling exponents C,p of the axis-perpendicular longitudinal velocity struc- 
ture functions, Sp = {\[v{r) — v{r + $y\ -^/^l^) ^ their extended self-similarity (ESS) 
first observed by Benzi et al. ( 1993[ ) in the non-rotating case is exploited. By this heuris- 
tic but generally accepted procedure exponents up to order p = 8 can be determined 
with sufficient precision in spite of the limited spatial resolution and overall duration 
of the simulations. The velocity field is resolved into a 2D component (fcy = 0) and 3D 
fiuctuations (fcy ^ 0) (see, e.g., Bourouiba fc Bartello|2007 1 . The relative exponents Cp/C2 
obtained via ESS from the 3D velocity coincide witn the Qp since the numerical results 
Ek±_ ~ ~ k'^^^'^'^^^ yield Cl^ = 1- For the 2D component = 2 is chosen consis- 
tently with the inertial-range scaling of the 2D energy spectrum ^ fcj^ (dotted line in 
figure [2|. 

The ESS-results shown in figure |4] and table |2] indicate a weakly intermittent structure 
of the 2D velocity component. The difference between exponents from non-rotating tur- 
bulence and those for f2 = 50 indicates that the reduction of intermittency for higher Q. 
is not a simple consequence of the removal of fluctuations with higher /cy . The 3D fluctu- 
ations exhibit an intermittent signature that is up to order p = 7 virtually independent 
of the rotation rate. The highest order scaling exponent is, however, suggesting slightly 
higher intermittency for = 5 than for = 50. 

A comparison with two-point statistics obtained with the total turbulent velocity fleld 
including all fcy contributions ( Miiller fc Thiele|2007 ) reveals that the 2D-component of 
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p 


Sp 


A2D,n=50 

Sp 


^nonintcrmitt 
Sp 


A3D,n=5 
Sp 


A3D,n=50 
Sp 


^modcl 
Sp 


1 


1.02 ±0.03 


1.0 ±0.01 


1 


0.52 ±0.003 


0.52 ±0.002 


0.53 


2 


2 


2 


2 


1 


1 


1 


3 


2.95 ±0.06 


2.99 ±0.02 


3 


1.43 ±0.006 


1.43 ±0.007 


1.41 


4 


3.88 ±0.1 


3.97 ±0.05 


4 


1.82 ±0.02 


1.80 ±0.02 


1.78 


5 


4.78 ±0.2 


4.94 ±0.08 


5 


2.15 ±0.03 


2.13 ±0.03 


2.11 


6 


5.67 ±0.3 


5.91 ±0.1 


6 


2.42 ±0.05 


2.42 ±0.05 


2.41 


7 


6.53 ±0.4 


6.86 ±0.1 


7 


2.65 ±0.07 


2.68 ±0.07 


2.68 


8 


7.38 ±0.5 


7.81 ±0.2 


8 


2.84 ±0.1 


2.91 ±0.09 


2.94 



Table 2. Structure function scaling exponents for 2D and 3D components of the velocity field 
with = 5 and — 50. The values are obtained by ESS assuming (^1° = 1 a^nd — 2. 



velocity tends to govern structure function scaling for = 50. In contrast, the Cp taken 
with = 5 are rather close to the values obtained with the fc|| velocity. The domi- 
nance of the 2D fluctuations in structure function scaling at high il is plausible because 
of the reduced activity of axis-perpendicular nonlinear interactions that is visible, e.g., 
in the corresponding energy flux ( .Miiller fc Thiele^^2007, ). The accompanying reduction 
of intermittency in the 2D flow component with higher rotation rate points towards 
stronger dynamical influence of inertial waves since weak wave-turbulence is generally 
non-intermittent. This effect is stronger in the 2D velocity than in the 3D component 
as the nonlinear energy flux is strongly depleted in the axis-parallel direction at high 
rotation rate compared to the axis-perpendicular flux (ibid.). 

The 3D structure function scaling is well represented by an intermittency model based 
on the She-Leveque (SL) ansatz ( She fc Leveque|1994[ ). The SL-model is quite successful 
in describing intermittent scalings of two-point statistics in various non-rotating tur- 
bulent flows. It can be written in a general form exhibiting two parameters which can 
be determined by physic al considerations ( Grauer et a?.||1994 Politano fc Pouquet||1995 
Miiller fc Biskamp|2000| : Cp = p/(35) + Co(l - (1 - 2/(3C^))^9). While Cq stands for the 



co-dimension of the most singular dissipative structures, g is connected to the basic non- 
intermittent scaling exponent, vi ~ £^/^ . In the present three-dimensional simulations the 
most-singular structures are quasi-one-dimensional vortex fllaments, i.e. Cq — 2, while 
relation (5.1) yields g=2. This results in the intermittency model 

Cp=p/6±2(l-(2/3)P/2) (3.1) 

which agrees well with the numerical data as shown in figure |4] (see also table [2]) . 

To conclude the analysis of the structure of rotating turbulence we present probability 
density functions (PDFs) P{Sv^) of the axis-perpendicular longitudinal velocity incre- 
ments Sv^ — [v{r) — v{r + ^)] • ^/^. Characteristic one-time PDFs for the 3D velocity 
field are shown in figure |5] For both rotation rates the PDFs shift from Gaussian form at 
large scales to functions with stretched exponential tails at small scales. This well-known 
behaviour indicates that both flows are intermittent in agreement with the structure 
function scaling. The tendency towards nonintermittent quasi-two-dimensional flow for 
increasing rotation rate is visible in the slight contraction of the normalized PDFs that 
is quantified by decreasing fiatness, 34/(82)'^, taking on values of 9.4 (fi = 5) and of 7.8 
(f] = 50) at small scales ^. At large scales the fiatness is, for both values of fi, close to 3, a 
value characteristic for a Gaussian distribution. The skewness 83/82^^ is also consistent 
with Gaussian statistics at large scales where it is ranging around zero. For = 50 this 
behaviour is observed at all scales. For Q = 5 the skewness drops to a level close to —0.5 
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Figure 4. Axis-perpendicular structure function scaling exponents C,p of velocity fluctuations 
with fc|| = (left) and fcy 7^ (right) up to order 8 for moderate, strong, and no r otation. Solid: 
n = 5, dotted: Q. — 50, dash-3-dot: 0,-0, dashed: intermittency model (3.1 1, dash-dotted: 
non-intermittent scaling, (p = p/2 {Q — 5) and (p = p {^1 = 50). 



at the small-scale end, a value known from non-rotating flow. Time-averaged flatness 
and skewness of the 3D velocity component are shown in figure |7] These findings suggest 



the existence of a self-similar limiting case for very high rotation rates (cf. Baroud et at. 



20021. 



Exemplary one-time PDFs for the 2D component of the velocity field are shown in 
figure |6] the associated time-averaged skewness and flatness functions are displayed in 
figure [Sj In agreement with the weakly intermittent scaling of the corresponding structure 
functions, the PDFs at all spatial scales are close to Gaussian with the flatness increasing 
for growing Q. from 2.5 (fi = 5) to 3.1 (fi = 50) at large scales and decreasing from 3.5 
(n = 5) to 3.2 (ri = 50) at small scales. With increasing rotation rate the skewness also 
changes from —0.1 (small scales) and about (large scales) for = 5 to positive values 
close to at all scales for Q. = 50. 



4. Turbulent energy decay under rotation 

The influence of rotation on the decay of turbulent energy still remains a controversial 
and yet unresolved problem. To help clarify this issue the time evolution of kinetic energy 
in freely decaying turbulence has been investigated in simulation III for eight different 
rotation rates ranging from = up to 51 = 5. 

Before describing the decaying runs, the macroscopic effects of the sudden onset of 
rotation on the kinetic energy in simulations I and II are reported. In these cases, E 
displays a sharp drop of about 20% (I) and 13% (II) with a subsequent remount that 
levels off close to the previous state. The dissipation rate e also drops but does not 
increase again, remaining at e ~ 0.05 in both simulations. This transient behaviour can 
be understood by the rotation-induced depletion of the spectral energy transfer that has 
been introduced earlier in combination with the applied forcing mechanism. For a quasi- 
stationary state of the turbulence the time-evolution of the total energy is characterized 
by £' = £/— e = 0, where £/ is the energy flux injected through the forcing. It is important 
to note that the forcing-mechanism used in simulations I and II does not provide a 
constant energy injection but rather an that nonlinearly adjusts to the "demand" of 
the fluctuations with k > kf and is thus coupled to the cascade dynamics. When rotation 
is switched on, the cascade gets depleted and consequently Sf diminishes. On the other 
hand, as less energy reaches the dissipative scales, e also decreases. These two effects lead 



Rotating homogeneous turbulence 



11 




Figure 5. PDFs of 3D velocity differences (fcy 7^ 0) taken at different spatial distances (increas- 
ing from left to center in each plot) for 57 = 5 (left) and f2 = 50 (right). The width is rescaled by 
the velocity variance, and all curves are shifted to a common maximum at 1 to allow for simpler 
comparison with the shape of a Gaussian (dash-dotted curve). 
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Figure 6. PDFs of 2D velocity differences (fcy = 0) taken at different spatial distances for n — 5 
(left) and — 50 (right). The width is rescaled and all curves are shifted as in figure |5] The 
dash-dotted curve represents a Gaussian. 



to a transient evolution of the energy E until a new equilibrium of forcing, cascade and 
dissipation has been reached. The observed behavior does not differ qualitatively when 
the rotation is ramped up (as has been checked by test computations). A similar drop of 
the initial energy has been reported by Yeung & Zhou ( 1998 1. There, however, the energy 
did not remain constant afterwards but increased until the end of the simulation. This 
discrepancy is probably due to their stochastic driving mechanism leading to a pile-up 
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Figure 7. Time-averaged flatness Si{Ci / {S2{£,)f (left) and skewness S-i{i) / {S2{C)f^'^ (right) 
of the PDFs corresponding to the 3D velocity fluctuations (figure [5| for for SI = 5 (solid line), 
f2 = 50 (dashed line), and Q. = Q (dash-3-dotted line). The dash-dotted lines display the respec- 
tive values for a Gaussian. 




Figure 8. Time-averaged flatness SiiC) / {S2{£,)f' (left) and skewness S-j,{i) / {S2{i)f^ (right) 
of the PDFs corresponding to the 2D velocity fluctuations (figure™ for for SI = 5 (solid line), 
f2 = 50 (dashed line), and $1 = (dash-3-dotted line). The dash-dotted lines display the respective 
values for a Gaussian. 



of energy at the largest scales. Similarly, the Gaussian forcing spectrum that ISmith & 



WalefFe ( 1999 ) applied at intermediate scales led to a growth of the total energy. 



In the decaying case III rotation is switched on after 1 LET of free evolution, i.e. t = 3. 
The initial random state with a Gaussian energy spectrum is similar to that of the driven 
systems I and II and the same for all values of fi. The temporal evolution of E{t) shown 
in figure |9] clearly demonstrates that the energy in rotating turbulence decays slower than 
in the non-rotating flow. This effect, a consequence of the attenuation of the spectral flux. 



has also been observed experimentally (see Jacquin et a?. [1990 ) and in LES (see Squires 
et aZ.||1993| [Yang fc Domaradzki|[2004t 



Since the modes with fcy = play a special role in rotating turbulence, it is instructive 
to regard also the 2D kinetic energy i?2D(^) — J dk± E{k\\ — 0,k±,t). As shown in figure 



10 i?2D(i) tends to be almost constant for higher rotation rates indicating a separation 



in the energetic evolution of 2D and 3D flow components (see, e.g., Bourouiba & Bartello 



20071 that becomes significant for 51 > 1.5. This trend does, however, not lead to total 



decoupling of 2D and 3D velocity components even in the limit fi ^ oo (Babin et al. 



20001. Consequently and since the relative contribution of 2D energy compared to 3D 
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Figure 9. Decay of the total kinetic energy for f7 = (dashed Hne) and Q, = 0.25, 0.5, 1, 1.5, 
2.5, 3.5 and 5 (straight hnes from bottom to top). 
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Figure 10. Decay of the 2D kinetic energy E2o{t) = J dki_ i5(fc|| = 0, k±,t) for O = (dashed 
hne) and Q = 0.25, 0.5, 1, 1.5, 2.5, 3.5 and 5 (straight lines from bottom to top). 



energy in the present simulations is also small, the decay of total energy in the present 
simulations is dynamically governed by the evolution of 3D turbulent fluctuations. Hence 
for simplicity and to facilitate comparison with other works on energy decay in rotating 
turbulence, the analysis of decay properties presented in the following focuses on the 
(2D+3D) kinetic energy. 

For all rotation rates, the energy exhibits a period of approximate self-similar decay, 
E ^ t~°', see figure [o) In isotropic Navier-Stokes turbulence this property has been 

), and it is of considerable interest to investigate 
how it is affected by rotation. During the power-law decay the logarithmic slope decreases 



extensively studied (see e. g. Lesiexir|1997 
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Figure 11. Logarithmic derivative of the total kinetic energy for Q. = Q (dashed line) and 
SI = 0.25, 0.5, 1, 1.5, 2.5, 3.5 and 5 (straight lines from bottom to top). 



n 0.25 0.5 1.0 1.5 2.5 3.5 5.0 

t™,„ 11.0 16.4 10.0 11.0 9.0 7.4 7.4 6.7 
t^a:. 20.1 33.1 22.2 18.2 20.1 16.4 16.4 14.9 

Table 3. Time intervals used for linear fit in figure |9]. 



with increasing rotation rate relating an fi-dependent decay exponent a. Earlier LES and 
weak turbulence computations seem to indicate that a is independent of the rotation rate 
Q. with a ~ ao/2, where ao is the non-rotating value (see Squires et a?.|1993" Bellet et al. 



20061. Theoretical predictions of Canuto & Dubovikov (1997) and Squires et al. (19931 



based on dimensional and scaling analysis also suggest a scaling exponent that is constant 
for all in contrast to our numerical results. 

The actual values of an were obtained by determining the slopes of linear fits to the 
approximately self-similar regions in the logarithmic plot of the energy. The respective 
time intervals were chosen to match the quasi-constant regions in the logarithmic deriva- 
tive shown in figure [TT] and are listed in table [3] The different lengths and starting points 
of the fit intervals produce some error in the read-off, that can be estimated by the 
change of the exponents when determined for slightly shifted intervals. Finally, we ob- 
tain an ri-dependence of the decay exponent as shown in figure |12[ For f2 = we find 
OfQ = 1-54 ± 0.16 in rough agreement with experimental and theoretical decay studies 



of non-rotating Navier-Stokes turbulence by Kolmogorov (19411, Saffman (19671, and 



Mohamed & LaRue ( 1990 1 . It should be noted that the absolute values of the 3D en- 



ergy decay exponents are systematically increased by about 0.4 as compared to the total 
3D-I-2D energy decay while exhibiting qualitatively the same f2-dependence (within error 
margins). This is a consequence of the additional quasi-constant contribution of the 2D 
fluctuations in the 3D-I-2D energy. 

The dependency of the decay exponent on the rotation frequency can be explained 
phenomenologically based on the scaling properties of rotating turbulence as provided 
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Figure 12. Decay exponents a{Q) for self-similar energy decay in rotating turbulence. The 
dashed line corresponds to a{Q) = ^j*^"^^ with ao — a{Q = 0) and tq ~ 0, 9. 



above. First, E ~ 1'°'^"^ and E 



e yields dimensionally 

E 

a — ^ £ 

To 



(4.1) 



with To being the large-eddy turnover time. Evaluating (5.11 at Lq, where most of the 
total energy E resides, gives 



Lin 



Here, we have used the relation Lq 

follow a{il) ~ l/{TQil). To include the correct asymptotics for 
end up with 



E2 tq. Inserting (4.2) in (4.1), a(ri) is found to 

0, i. e. a 



OLn 



or in a shorter notation, 



ao 



1 



ao, we 



(4.3) 



(4.4) 



Figure [12] indicates that this simple model qualitatively reproduces the general features 
of the fi-dependence of the scaling exponent. This fi" ^-dependence suggests a complete 
depletion of nonlinear decay in the limit U, ~^ oo. Such behaviour would also be in accord 
with a high level of two-dimensionalization, i.e. E2D ^ £-30 • We note, however, that the 
phenomenology presented here will probably cease to be valid in this asymptotic limit 
where weak turbulence theory should apply. 



5. Conclusion 

In summary, high-resolution direct numerical simulations of large-scale-driven and de- 
caying rotating hydrodynamic turbulence have been conducted for slow and rapid rota- 
tion and analyzed anisotropically. With increasing rotation rate, all considered systems 
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show a trend toward dynamical two-dimensionalization in planes perpendicular to the 
rotation axis, strong attenuation of the nonlinear spectral energy flux along $7, and con- 
centration of energy around the plane fcy = 0. The modification of the energy cascade 
by rotation leads to an energy spectrum that scales as ~ k~[^ in the incrtial range with 
no discernible scaling in /cy . The perpendicular scaling results from the integration over 
spectral contributions for all fey which exhibit A:~'^-behaviour near fcy = for low ro- 
tation rate. Also a modified approximate self-similar decay of the total kinetic energy 
E ~ is observed with a ^ 17"^. Separating the velocity field into 2D {k\\ = 0) and 

3D (fc|| ^ 0) components reveals that the 2D energy tends to become time-independent 
as the rotation rate increases. Energy decay and spectral scaling can be reproduced by 
simple phenomenological models which also yield a non-intermittent scaling relation for 
fluctuations perpendicular to the rotation axis, 

~ (efi)V4^V2 , (5.1) 

with ^ being an axis-perpendicular scale. The structure function scaling and the PDFs 
show weak intermittcncy of the 2D velocity component while the intermittent signature 
of the 3D fluctuations is well represented by a model based on the She-Leveque ansatz. 
The intermittcncy of the 3D velocity component is reflected in the respective PDFs by 
the well-known change from Gaussian at large scales to leptocurtic at small scales. The 
two-point statistics show an overall weak dependence on the rotation rate. 

The authors thank Stephan Kiimmel, Walter Zimmermann and Lorenz Kramer for 
their support. WCM gratefully acknowledges discussions with J. Leorat, A. Mahalov, F. 
Moisy and J. Seiwert. 
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